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Abstract. We show that the space of Euclid's parameters for Pythagorean 
triples is endowed with a natural symplectic structure and that it emerges as 
a spinor space of the Clifford algebra R21, whose minimal version may be con- 
ceptualized as a 4-dimensional real algebra of "kwaternions." We observe that 
this makes Euclid's parameterization the earliest appearance of the concept of 
spinors. We present an analogue of the "magic correspondence" for the spinor 
representation of Minkowski space and show how the Hall matrices fit into 
the scheme. The latter obtain an interesting and perhaps unexpected geomet- 
ric meaning as certain symmetries of an Apollonian gasket. An extension to 
more variables is proposed and explicit formulae for generating all Pythagorean 
quadruples, hexads, and decuples are provided. 
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1. Introduction 

In the common perception, Clifford algebras belong to modern mathematics. It 
may then be a point of surprise that the first appearance of spinors in mathematical 
literature is already more than two thousand years old! Euclid's parameters m and 
n of the Pythagorean triples, first described by Euclid in his Elements, deserve the 
name of Pythagorean spinors — as will be shown. 

The existence of Pythagorean triples, that is triples of natural numbers (a, 6, c) 
satisfying 



has been known for thousands of year^J Euclid (ca 300 BCe) provided a formula 
for finding Pythagorean triples from any two positive integers m and n, m > n, 



1 The cuneiform tablet known as Plimpton 322 from Mesopotamia enlists 15 Pythagorean triples 
and is dated for almost 2000 BCE. The second pyramid of Giza is based on the 3-4-5 triangle 
quite perfectly and was build before 2500 BCE. It has also been argued that many megalithic 
constructions include Pythagorean triples |17| . 





2 



Jerzy Kocik 



namely: 



2 2 

a — m — n 



b = 2 ran (1.2) 

2 i 2 

c — m + n 

(Lemma 1 of Book X). It easy to check that a, 6 and c so defined automatically 
satisfy Eq. (1.1), i.e., they form an integer right triangle. 

A Pythagorean triple is called primitive if (a, 6, c) are mutually prime, that is 
gcd (a, 6, c) = 1. Every Pythagorean triple is a multiple of some primitive triples. 
The primitive triples are in one-to-one correspondence with relatively prime pairs 
(m, n), gcd (m, n) = 1, m > n, such that exactly one of (m, n) is even [15} [T6]. 

An indication that the pair (m, n) forms a spinor description of Pythagorean 



triples comes from the well-known fact that (1.2) may be viewed as the square 
of an integer complex number. After recalling this in the next section, we build a 
more profound analysis based on a 1:2 correspondence of the integer subgroups of 
0(2, 1) and 51^(2, R), from which the latter may be viewed as the corresponding 
pin group of the former. For that purpose we shall also introduce the concept of 
pseudo-quaternions ( "kwaternions" ) — representing the minimal Clifford algebra 
for the pseudo-Euclidean space M 2 ' 1 . 

The most surprising result concerns an Apollonian gasket where all the objects 
mentioned acquire a geometric interpretation. 



2. Euclid's Labels as Complex Numbers 

Squaring an integer complex number z = m+ni will result in an integer complex 
number 

z 2 = (m + ni) 2 — (m 2 — n 2 ) + 2mn z, (2.1) 
the norm of which turns out to be also an integer: 

\z 2 \ = m 2 + n 2 . (2.2) 

This auspicious property gives a method of producing Pythagorean triangles: just 
square any integer complex number and draw the result. For instance z = 2 + i will 



produce 3-4-5 triangle, the so-called Egyptian triangle (Figure |2p). Equation 2.1 
is equivalent to Euclid's formula for the parameterization of Pythagorean triples 
[T6] . We shall however allow z to be any integer complex number, and therefore 
admit triangles with negative legs (but not hypotenuses). The squaring map 

sq: C C : z -> z 2 

(exclude zero) has a certain redundancy — both z and - z give the same Pythagorean 
triple. This double degeneracy has the obvious explanation: squaring a complex 
number has a geometric interpretation of doubling the angle. Therefore one turn 
of the parameter vector z = m + ni around the origin makes the Pythagorean 
vector z 2 go twice around the origin. We know the analogue of such a situation 
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(m+n/)^ 




(b) 



z = 2 + i 




z 2 =3 + 4/ 



Figure 1. (a) Euclid's parameterization, (b) The Egyptian trian- 
gle is the square of z = 2 + i. 

from quantum physics. Fermions need to be turned in the visual space twice be- 
fore they return to the original quantum state (a single rotation changes the phase 
of the "wave function" by 180°). Mathematically this corresponds to the group 
homomorphism of a double cover 

SU{2) ^ 50(3) 

which is effectively exploited in theoretical physics [2]. The rotation group SO(3) 
is usually represented as the group of special orthogonal 3x3 real matrices, and 
the unitary group SU(2) is typically realized as the group of special unitary 2x2 
complex matrices. The first acts on R 3 — identified with the visual physical space, 
and the latter acts on C 2 — interpreted as the spin representation of states. One 
needs to rotate a visual vector twice to achieve a single rotation in the spinor space 
C 2 .This double de generacy of rotations can be observed on the quantum level pQ, 



Euclid's 
parameter 
space 




Space of 
Pythagorean 
vectors 



Figure 2. Double degeneracy of the Euclid's parameterization. 

Remark 2.1. The group of rotations SO{3) has the topology of a 3-dimensional 
projective space, and as such has first fundamental group isomorphic to Z2. Thus 
it admits loops that cannot be contracted to a point, but a double of a loop 
(going twice around a topological "hole") becomes contractible. Based on this 
fact, P.A.M. Dirac designed a way to visualize the property as the so-called "belt 
trick" [ME]. 
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This analogy between the relation of the Pythagorean triples to the Euclid's 
parameters and the rotations to spinors has a deeper level — explored in the 
following sections, in which we shall show that (2.1) is a shadow of the double 
covering homomorphism 5X^(2) — >> 0(2, 1). It shall become clear that calling the 
Euclid parameterization (m, n) a Pythagorean spinor is legitimate. 

Before we go on, note another interesting redundancy of the Euclidean scheme. 
For a complex number z = m + ni define 

z d = (m + n) + (m — ri)i . 

The square of this number brings 

(z d f = Amn + 2(m 2 - n 2 ) i (2.3) 

with the norm \z d \ 2 = 2(m 2 + n 2 ). Thus both numbers z and z d give - up to scale 
- the same Pythagorean triangle but with the legs interchanged. For instance both 
z — 2 + i and z d = 3 + i give Pythagorean triangles 3-4-5 and 8-6-10, respectively 
(both similar to the Egyptian triangle). In Euclid's recipe (1.2), only one of the 
two triangles is listed. This map has "near-duality" property: 

(z d ) d = 2z. (2.4) 



3. Pseudo-quaternions and Matrices 

A. Pseudo-quaternions. Let us introduce here an algebra very similar to the algebra 
of quaternions. 

Definition. Pseudo-quaternions K (kwaternions) are numbers of type 

g = a + 6i + cj + dk (3.1) 

where a, 6, c, d E K and where i, j, k are independent "imaginary units". Addition 
in K is defined the usual way. Multiplication is determined by the following rules 
for the "imaginary units" : 

i 2 = i y = k 

j 2 = l jk = -i (3.2) 

k 2 = -1 ki = -j 

plus the anticommutation rules for any pair of distinct imaginary units: ij = — ji, 
jk = — kj, and ki = — ik. [The rules are easy to remember: the minus sign appears 
only when k is involved in the product]. 

The pseudo-quaternions form an associative, non-commutative, real 4-dimensional 
algebra with a unit. Define the conjugation of q as 

q = a — hi — cj — dk (3.3) 

and the squared norm 

\q\ 2 = qq = a 2 - b 2 - c 2 + d 2 (3.4) 
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The norm is not positive definite, yet kwaternions "almost" form a division algebra. 
Namely the inverse 

q- l =q/\q\ 2 (3.5) 

is well-defined except for a subset of measure zero — the "cone" a 2 — b 2 — c 2 +d 2 = 0. 

Other properties like \ab\ 2 = |6a| 2 , |a6| 2 = |a| 2 |6| 2 , etc., are easy to prove. The 
set G — {q E K | \q\ 2 = 1} forms a group, and G Q = SX(2,R). Topologically, the 
group (subset of kwaternions) is a product of the circle and a plane, G = S 1 x R 2 
. The bigger set G — {q E K | \q\ 2 = ±1} is isomorphic to S ± L(2 1 i?), the modular 
group of matrices with determinant equal ±1. 

The regular quaternions describe rotations of R 3 . Pseudo-quaternions describe 
Lorentz transformations of Minkowski space R 2 ' 1 , 0(2,1). Indeed, note that the 
map 

v —> v — qvq~ x (3-6) 

preserves the norm, \qvq~ 1 ] 2 = \v\ 2 . Thus the idea is quite similar to regular 
quaternions. Represent the vectors of space-time R 2,1 by the imaginary part 

v = xi + yj + tk (3.7) 

Its norm squared is that of Minkowski space R 2 ' 1 , |v| 2 = —x 2 — y 2j rz 2 . Define 
a transformation kwaternion as an element q E G. Special cases are: 

rotation in x?/-plane: q = cos 0/2 + k sin 0/2 

hyperbolic rotations q = cosh 0/2 + isinh0/2 (boost in direction of y-axis) 
(boosts) : q = cosh + j sinh 0/2 (boost in direction of x-axis). 

(3.8) 

It is easy to see that the norm of v is preserved transformation (3.7), as well 
as its lack of real part. 

B. Pseudo-quaternions as a Clifford algebra. The algebra of pseudo-quaternions 
is an example of a Clifford algebra. Namely, it is the Clifford algebra of a pseudo- 
Euclidean space R 2,1 , a 3-dimensional Minkowski space with quadratic form 

g = -x 2 - y 2 + z 2 (3.9) 

The corresponding scalar product of two vectors will be denoted as g(v, w) = 
(v, w). Let {fi, f2, fa} be an orthonormal basis of R 2 ' 1 such that: 

(fi,fi) = -i 

(f 2 ,f 2 ) = -i , N 

(f3,f3> = 1 (310) 

and (fj, fj) = for any i ^ j . 

Let us build the Clifford algebra [13] of this space. We shall assume the stan- 
dard sign convention: 

vw + wv = — 2(v, w) 
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(Clifford products on the left side, pseudo-Euclidean scalar product on the right). 
The basis of the universal Clifford algebra R2,i of the space R 2,1 consists of eight 
elements 

R 2 ,i =span{l, fi, f 2 , f3, fif2, f2f3, f3fi, f&k} 
The basic elements satisfy these relations (in the sense of the algebra products (cf. 
B) : 

1 



fl 2 
ff 

f 2 

and fifj 



1 



(3.H) 



-ijfi for any i ^ j . 

Relations between the other elements of the basis of the Clifford algebra are in- 
duced from these via the associativity of the algebra product. 

C. Matrix representation of pseudo-quaternions. Note that the following four 



matrices satisfy the above relations (3.11): 

























" 1 


" 




" 


1 " 




"-10" 




"01" 


cr = 





1 


£71 = 


1 





02 = 


1 


0"3 = 


-1 



(3.12) 



4 



4 



4 



-<7 . 



(3.13) 



(3.14) 



Indeed, the squares are, respectively: 

2 

The products are 

G\G2 — +^"3 CT2^"3 — — —&2 

and o§Oi — (J{ for each i. Thus these matrices play a similar role as the Pauli 
matrices in the case of spin description or the regular quaternions. Note however 
the difference: the above construct is over real numbers while the Pauli matrices 
are over complex numbers. 

The other elements are represented by matrices as follows: 

1 =cr fi f2 cr3 fa f 1 —&2 fihk — ^0 • 

As in the case of quaternions, the universal Clifford algebra for R 2,1 is 8-dimensional, 
yet the minimal Clifford algebra is of dimension 4 and is isomorphic to the algebra 
of kwaternions K. 

Remark 3.1. In algebra, pseudo-quaternions are also known as split quaternions, 
especially in the context of the CayleyDickson construction [3j. Other names in- 
clude para-quaternions, coquaternions and antiquaternions. 
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4. Minkowski Space of Triangles and Pythagorean Spinors 

Now we shall explore the geometry of Euclid's map for Pythagorean triples. As a 
map from a 2-dimensional symplectic space to a 3-dimensional Minkowski space, 
we have: 

^:(M 2 ,a;)^(R 2 >\G9 

(m, n) — >► (m 2 — n 2 , 2mn, m 2 + n 2 ) 



The map (2.1) of squaring complex numbers is just its truncated version. We shall 
discuss symmetries of the natural structures of both spaces. Obviously, we are 
mostly interested in the discrete subsets 1? and Z 3 of those spaces. 

First, introduce the space of triangles (x, y, z) as a real 3-dimensional Minkowski 
space M 2 ' 1 with a quadratic form 



Q = -x 2 



2 , 2 

y +z 



(4.2) 



("space-time" with the hypotenuse as the "time"), and with metric given by a 
3x3 matrix G — diag(l, 1,-1). The right triangles are represented by "light-like" 
(null) vectors, and the Pythagorean triples by the integer null vectors in the light 
cone. Clearly, not all vectors correspond to real triangles, and the legs may assume 
both positive and negative values. We have immediately 

Proposition 4.1. The group of integer orthogonal matrices 0(2, 1;Z) C 0(2, 1;M) 
(Lorentz transformations) permutes the set of Pythagorean triangles. 

On the other hand we have the two-dimensional space of Euclid's parameters 
E = R 2 . Its elements will be called Pythagorean spinors. Occasionally we shall use 



the isomorphism 



P 2 ~ 



and identify [m, n] = m + ni. 



The space of Euclid's parameters will be equipped with an inner product: 

Definition 4.2. For two vectors u = [m,n] T and w = [ra / ,n / ] T , the value of the 
symplectic form uj is defined as 

cj(u, w) — run — nm! . (4.3) 

Conjugation A* of a matrix A representing an endomorphism in E is the adjugate 
matrix, namely 



if A 



b 
d 



then A* 



d 



a 



(4.4) 



Conjugation of vectors in E is a map into the dual space, expressed in terms of 
matrices as 

r rn 



n 



[- 



n m 



(4.5) 



Now, the symplectic product may be performed via matrix multiplication: u(u, w) = 
u*w. The map defined by ( |4.5[ ) is the symplectic conjugation of the spinor. Also, 
note that AA* = A* A = det(A)J. 

Remark 4.3. In the complex representation, the symplectic product is uj{u,w) — 
7}(uw — wu). 
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(0 


(AB)* 


= B*A* 


(it) 


(A)** 


= A 


(in) 


(An)* 


= u*A* 



Proposition 4.4. For any two matrices A and B and vector u in the spinor space 
we have: 



(4.6) 



Proposition 4.5. The group that preserves the symplectic structure (up to a sign) 
is the modular group SL ± {2^ T L) C 5X^(2, R) understood here as the group of '2x2 
integer (respectively, real) matrices with determinant equal ±1. 

Proof. Preservation of the symplectic form is equivalent to matrix property A* A = 
±1. Indeed: 

u{Au,Aw) (Au)*(Aw) u(A*A)w det(A)u*w ±u(u,w) 

for any u,w, since A* A = ± det (A) 7", and by assumption det (A) = ±1. Let us 
relate the two spaces. Given any endomorphism M of the space of triangles R 2 ' 1 , 
we shall call an endomorphism M its spinor representation, if 

M(<p(u)) = <p(Mu) . (4.7) 

Now we shall try to understand the geometry of the spin (Euclid's) representation 
of the Pythagorean triples in terms of the kwaternions defined in the previous 
sections. Recall the matrices representing the algebra: 

^0 

By analogy to the spinor description of Minkowski space explored in theoretical 
physics, we shall build a "magic correspondence" for Pythagorean triples and their 
spinor description. First, we shall map the vectors of the space of triangles, R 2 ' 1 
into the traceless 2x2 real matrices, M® 2 • The map will be denoted by tilde 



" 1 


" 







1 ' 




'-10' 




1 " 





1 


(71 = 


1 





C 2 = 


1 


0"3 


-1 



p2,l 



M®2 and denned: 



v (x,y,z) 





1 


"01' 


+ y 


i 


" -1 ' 


1 


1 " 


— > V = x ■ 


2 


1 


' 2 


1 


+ Z -2 


-1 


1 




-y x + z 










~ 2 


X 


— z 


y 











(4.8) 



□ 

Proposition 4.6. The matrix representation of the Minkowski space of triangles 
R 2 ' 1 realizes the original Minkowski norm via the determinant: 

1 1 v| | = 4 det v — —x 2 — y 2 + z 2 . 

The scalar product may be realized by traces, namely: 

v • w = —2 Tr v w 

and particular vector coefficients may be read from the matrix by: 

v l — — det (Tj • Tr w ? - 
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The technique of Clifford algebras allows one to represent the orthogonal 
group by the corresponding pin group. Since an orthonormal transformation may 
be composed from orthogonal reflections in hyperplanes, one finds a realization of 
the action of the Lorenz group on the Minkowski space of triangles via conjugation 
by matrices of the spin group 5X^(2, K); in particular, for any orthogonal matrix 
A E 0(2, 1; Z), the action v' = Av corresponds to 



; ' = AvA* 



(4.9) 



that is, the following diagram commutes: 

V = \ Y1 vi(J i 

conj A 



v 

A 

Aw 



And now the reward: since the Pythagorean triples lie on the "light cone" of the 
Minkowski space, we may construct them from spinors in a manner analogous to 
the standard geometry of spinors for relativity theory. But here we reconstruct 
Euclid's parameterization of the triples. Recall that 



X 




m 2 - n 2 


y 




2mn 


z 




m 2 + n 2 



Thus we have: 

Theorem 4.7. The spin representation of Pythagorean triples splits into a tensor 
(Kronecker) product: 

,2 1 [ <m 1 

[ 



—y x + z 
x — z y 



-mn rn 

.2 



-n 



mn 





m 




n 



-n rn 



(4.10) 



Hence we obtain yet another aspect of the Euclid's formula, namely a tensor 

(4.11) 



version of the (1.2) and (4.1): 



<p(u) 



u($ u 



Note that [— n,m] in (4.10) is the symplectic conjugation of the spinor [m,n] T . 
Now, due to the above Theorem, the adjoint action splits as follows 



rn 
n 



[-n m) A* = ( A 



m 




m 


n 




n 



(4.12) 



And the conclusion to the story: The spin representation emerges as "half of the 
above representation: 



rn 
n 



m 
n' 



= A 



m 
n 



(4.13) 
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Magic Correspondence 




Minkowski space M 2 ' 1 


Traceless 2x2 matrices 


Main object 

Norm 

Action 

Minkowski scalar product 
The ith coefficient 


v = (x,y,z) 
|| v || = -x 2 -y 2 + z 2 

v' = Av 

(AeO(2,l;Z)) 

v • w = v T Gw 
= v • 


||v|| = 4detv 

(AeSL^Z)) 
v • w = — 2 Trvw 
v % = — det <7i • Tr wi 



Table 1 . Magic correspondence for Pythagorean triples and their 
spinor description 



Remark 4.8 (on d-duality). The "duality" (2.4) is also represented in spin language, 
namely, the matrix D of exchange of x with y and the corresponding 2-by-2 spin 
matrix, D, are: 



D 



10 

1 
1 



D 



1 1 
1 -1 



(4.14) 



Indeed, we have: 



" 1 1 




rn 




m + n 


1 -1 




n 




m — n 



The d-duality may be expressed in the form of a commuting diagram 







X 




rn 










> 


y 




n 










z 





rn — n 
2mn 
m? + n 2 



D 



D 



m + n 
m — n 





y 




Amn 


+ 2 


X 




2m 2 - 2n 2 




z 




2m 2 + 2n 2 
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5. Hall Matrices and Their Spinor Representation 

It is known that all primitive Pythagorean triples can be generated by the following 
three Hall matrices [I] 



(5.1) 





' 1 


2 


2 " 




' 1 


-2 


2 " 




' -1 


2 


2 " 


u = 


2 


1 


2 


L = 


2 


-1 


2 


R = 


-2 


1 


2 




2 


2 


3 




2 


-2 


3 




-2 


2 


3 



by acting on the initial vector v 
a few examples: 



[3,4, 5] T ("Egyptian vector"). For clarification, 



Lv 



Uw 



-2 2 
-1 2 
-2 3 

12 2 
2 12 
2 2 3 





" 3 " 




5 




4 




12 




5 




13 




" 3 " 




" 21 " 




4 




20 




5 




29 



Rv 



" -1 2 2 " 




" 3 " 




" 15 " 


-2 1 2 




4 




8 


-2 2 3 




5 




17 



URL 2 £/v 



3115 
3348 
4573 



Let us state it formally: 



Theorem 5.1 (Hall). The set of primitive Pythagorean triples is in one-to-one 
correspondence with the algebra of words over alphabet {U,L,R}. 

The original argument of Hall was algebraic (see j4j for a proof). But here we 
shall reinterpret this intriguing result in terms of geometry. Hall matrices may be 
understood in the context of our previous section and augmented with the spinor 
description. Let us start with this: 

Proposition 5.2. The Hall matrices and their products are elements of the Lorentz 
group 0(2, 1; Z). In particular, for any v = [x, y, z] T , 

j(v,v) = 0^ #(Mv,Mv) = 0, 

and therefore they permute Pythagorean triples. 

Proof. Elementary. Recall that the matrix of the pseudo-Euclidean metric is G = 
diag( 1,1,-1). One readily checks that Hall matrices preserve the quadratic form, 
i.e., that X T GX = G for X = U, L, R. In particular we have L, R E 50(2, 1), as 
det L = 1, det R = 1. Since det U = —1, U contains a reflection. □ 

Hall's theorem thus says that the set of primitive Pythagorean triples coincides 
with the orbit through [3,4, 5] T of the action of the semigroup generated by the 
Hall matrices, a subset of the Lorenz group of the Minkowski space of triangles: 

gen{i?,L,£/} C 0(2,1;Z). 

Thus the results of the previous section apply in particular to the Hall semi- 
group. In particular: 
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Theorem 5.3. The spin representation of the Hall matrices are 



u =[i o\ L =[i ~o J R =[o i\ ■ ^ 

In particular, they are members of the modular group 51/^(2, Z) C 5X^(2, R), that 
is they preserve the symplectic form up to a sign. 

Proof The transformations in the space of Euclid's parameters that correspond to 
the Hall matrices may be easily found with elementary algebra. By Proposition |45j 
we need to check M*M = ±7. Simple calculations show: 



t/*C7 = —J L*L = 7 R*R = I. 
That is, the determinants are det U = — 1, det L = 1, and det R 



1. 



□ 



Clearly, the "magic correspondence" outlined in the previous section holds as 
well for the Hall matrices. In particular, the spin version of matrices act directly 



on Pythagorean spinors (4.9), and the action on the Pythagorean vectors may be 



obtained by the tensor product (4.13) 



A classification of the semigroups in the modular group 5L ± (2, Z) that gen- 
erate the set of primitive Pythagorean spinors as their orbits seems an interesting 
question. 

Remark on the structure of the Pythagorean semigroup and its spin version. The 

last unresolved question concerns the origin or the geometric meaning of the Hall 
matrices and their spin version. First, one may try to interpret Hall matrices in 
terms of the Minkowski space-time structure. They may easily be split into a boost, 
spatial rotation and reflection. Indeed, define in R 2,1 these three operators: 





3 









c s 


" 


H = 





1 





, T = 


—s c 







_V8 





3 







1 



where c 



^2/2 



(H = boost by "velocity" (>/3/8,0) , T = rotation by 45°). Then the Hall matrices 
are the following Lorentz transformations : 

U = T 2 (T~ 1 HT) = THT 
R = THTRi 
L = THTR 2 

where R\ — diag( — 1,1,1) and i?2 = diag(l,— 1,1) represent reflections. For in- 
stance, U represents a boost in the special direction (1,1) followed by space point- 
inversion. This path seems to lead to nowhere. Thus we may try the spinor version 
of the Hall matrices. One may easily see that the latter can be expressed as linear 
combinations of our pseudo-Pauli basis: 



^"0 — &2 — C"3 i? = CTq + CTi + (J3 U — (Jq + G\ — <J<i 
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But this naive association does not seem to explain anything. Quite surprisingly, 
insight may be found in the geometry of disk packing. This is the subject of the 
next section. 

6. Pythagorean Triangles, Apollonian Gasket and Poincare 
Disk 

Now we shift our attention to — at first sight rather exotic for our problem — 
Apollonian gaskets. 

Apollonian window. Apollonian gasket is the result of the following construction. 
Start with a unit circle, called in the following the boundary circle. Inscribe two 
circles so that all three are mutually tangent. Then inscribe a new circle in every 
enclosed triangular-shaped region (see Figure 3a). Continue ad infinitum. A spe- 
cial case, when the first two inscribed circles are half the radius of the boundary 
circle will be called the Apollonian window (shown in Figure^). The Apollonian 




Figure 3. An Appollonian gasket and the Apollonian window 

window has amazing geometric properties [TT] . [9]. One of them is the fact that 
the curvature of each circle is an integer. Also, as demonstrated in [8], the center 
of each circle has rational coordinates. It can be shown that the segment that joins 
the centers of any two adjacent circles forms the hypotenuse of a Pythagorean 
triangle, whose two legs are parallel/perpendicular to the main axes. More specif- 
ically, the sides of the triangle become integers when divided by the product r\r2 
of the radii of the adjacent circles (see Figure [I]). 

Let us call a subboundary circle any circle in the Apollonian window tangent 
to the boundary circle (shaded circles in Fig. [5]). Consider a pair of tangent circles 
of which one is the boundary — and the other a subboundary circle. If we prolong 
the hypotenuses of the associated triangles, we shall hit points on the boundary 
circle, namely the points of tangency. Due to this construction, the slope of each 
such line is rational. We will try to see how to permute these points. 

Among the many symmetries of the Apollonian window are inversions in the 
circles that go through the tangency points of any three mutually tangent circles. 
Such inversions permute the disks of the window, and in particular preserve their 
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Figure 4. Pythagorean triangles in the Apollonian window. Bold 
numbers represent the curvatures of the corresponding circles. 

tangencies. We shall look at the following three symmetries labeled A, B, C (see 
Figure [5)3) : 

A - reflection through the vertical axis 

B - reflection through the horizontal axis (6.1) 
C - inversion through the circle C (only a quarter of the circle is shown) 

Proposition 6.1. The three compositions of maps 

CA, CB, CBA (6.2) 

leave the set of subboundary circles of the first quadrant invariant. In particular, 
they permute points of tangency on the boundary circle in the first quadrant. 

Proof. A reflection in line A or B or their composition AB (reflection through 
the central point) carries any circle in the first quarter to one of the other three 
quarters. If you follow it with the inversion through C, the circle will return to the 
first quarter. Since tangency is preserved in these transformations, the proposition 
holds. □ 

The crucial observation is that both the lines A, B and circle C may be 
understood as "lines" in the Poincare geometry, if the circle is viewed as the 
Poincare disk. This will allow us to represent these operations by matrices using 
the well-known hyperbolic representation of the Poincare disk. 

Poincare disk. Recall the geometry of the hyperbolic Poincare disk. In the standard 
model, the set of points is that of a unit disk D in the Euclidean plane: 

D = {(x, y) e M 2 : x 2 + y 2 < 1} with dD = {(x, y) : x 2 + y 2 = 1}. 
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Figure 5. (a) Subboundary circles and the corresponding 
Pythagorean rays, (b) Three symmetries of the Apollonian window. 



Poincare lines are the circles that are orthogonal to dD (clearly, only the intersec- 
tion with D counts). This geometry - as is well known - may be induced from a 
hyperbolic linear space. Consider a three-dimensional Minkowski space M 2 ' 1 and a 
hyperboloid H: 

t 2 - x 2 - y 2 = 1 . (6.3) 

Stereographic projection tt onto the plane P E R 2 ' 1 defined by t = 0, with the 
vertex of projection at (—1,0,0), brings all points of the hyperboloid H onto D 
in a one-to-one manner. In particular, each Poincare line in D is an image of 
the intersection of a plane in R 2 ' 1 through the origin O with the hyperboloid i7, 
projected by tt onto D. We shall use this plane-line correspondence. Recall also 
that reflection in a plane P can be done with the use of a unit normal vector n: 

R n :v -> v' = v-2Mn (6.4) 

(n, n) 

where the orthogonality n_LP and the scalar product are in the sense of the pseudo- 
Euclidean structure of the Minkowski space 



Back to the Apollonian window. Consider the three symmetries of the Apollo- 
nian window (6.2). Each of them may be realized in terms of a reflection in a 
corresponding plane in the hyperbolic representation. 
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Proposition 6.2. The three symmetries (6.1) of the Poincare disc (coinciding with 



the Apollonian window) have the following matrix representations 



Symmetry A 



Symmetry B 



Symmetry C 



ni = [1,0, 0] 5 



n 2 = [0,1, 0] 5 



n 3 = [1,1,1] 



T 



R^ = 



i?2 



i2a 



-1 











1 











1 


1 





" 





-1 











1 


-1 




2 2 


-2 




1 2 


-2 




2 3 



(6.5) 



Proof. One can easily verify that each is unit and corresponds to the assigned 
symmetry. Here are direct calculations for finding the matrix corresponding to the 
third symmetry C corresponding to n = 113. Acting on the basis vectors and using 



(6.4) we get: 



i? 3 ei = ei - 2 



i? 3 e2 = e 2 — 2 



(ei,n) 
(n,n) 



(e2,n) 
(n,n) 





" 1 " 




-1 


" 1 " 




" -1 " 













- 2 • 




1 




-2 











1 




-2 




" " 




-1 


" 1 " 




" -2 " 




1 


- 2 • 




1 




-1 











1 




-2 



Rse 3 = e 3 - 2 



(ei,n) 
(n,n) 



" " 


-2- — 


" 1 " 




" 2 " 





1 




2 


1 


-1 


1 




3 



which indeed defines matrix i? 3 . The other two, R\ and R2 are self-explanatory. □ 

And now we have our surprising result: 
Theorem 6.3. The hyperbolic representation of the permutations (6.2) of the sub- 



boundary circles of the first quadrant correspond to the Hall matrices: 



L 





' 1 


-2 


2 


R3R1 = 


2 


-1 


2 




2 


-2 


3 




" -1 


2 


2 


R3R2 — 


-2 


1 


2 




-2 


2 


3 




" 1 


2 


2 


R3R1R2 — 


2 


1 


2 




2 


2 


3 



= R 



U. 



(6.6) 
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Figure 6. Spinors for some Pythagorean triangles in the Apol- 
lonian window, given as labels at the points of tangency of the 
corresponding pairs of circles. 



Corollary. Each primitive Pythagorean triangle is represented in the Apollonian 
window. 



This finally answers our question of Section [5] on the origin and structure 
of the matrices of spin representation of the Hall matrices. Following the above 
definitions, we get in Clifford algebra the following representations of the reflections 



that constitute the symmetries (6.2): 



(fi + f 2 + f 3 )fi = + f 2 fi + fsfi = "(73 + - °2 
(fi + f 2 + f 3 )f 2 = fif 2 + f| + hh = (Jo + V3 + t7i = 

(fl + f 2 + f 3 )(-f 3 ) = "fife - f 2 f 3 - f 3 2 = -*2 + ^1 + ^0 



2 -1 
1 

1 2 
1 



R 



(6.7) 



2 1 
1 



U 



(cf. ( |5.2D ). Note that the duality map ( |4.13| ) that exchanges x and y in the 
Pythagorean triangles and which has a spinor representation D — o\ — a 2 fits 
the picture, too, as the vector n = [1,-1, 0] represents (scaled) reflection in the 
plane containing the x-y diagonal. 

We conclude our excursion into the spin structure of the Euclid parameteri- 
zation of Pythagorean triples with Figure [6] that shows spinors corresponding to 
some Pythagorean triangles in the Apollonian window. 
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7. Conclusions and Remarks 

We have seen that the space of Euclid's parameters for Pythagorean triples is 
endowed with a natural symplectic structure and should be viewed as a spinor 
space for the Clifford algebra ^2,1? built over 3-dimensional Minkowski space, 
whose integer light-like vectors represent the Pythagorean triples (see Figure [7]). 
The minimal algebra for R 2,1 is four-dimensional and may be conceptualized as 
"pseudo-quaternions" K and represented by 2 x 2 matrices. In this context the 
Pythagorean triples may be represented as traceless matrices, and Euclid's pa- 
rameterization map as a tensor product of spinors. This set-up allows us to build 
the spinor version of the Hall matrices. The Hall matrices acquire a geometric in- 
terpretation in a rather exotic context of the geometry of the Apollonian window. 

Euclid's discovery of the parameterization of Pythagorean triples may be 
viewed then as the first recorded use of a spinor space. The spinor structure of 
the Apollonian window is another interesting subject that will be studied further 
elsewhere. 

The method may be generalized to other dimensions, as indicated here: 

Method A. Define a Pythagorean (fc, /)-tuple as a system of (k + I) integers that 
satisfy 

a 2 + a\ + a 2 + . . . + a\ = b\ + b\ + b\ + . . . + b] (7.1) 
In order to obtain a parameterization of Pythagorean (fc, /)-tuples do the follow- 
ing: Start with pseudo-Euclidean space M. k >\ build a representation of the Clifford 
algebra Rfc,z . Then split the matrix that represents the isotropic vectors of M, k > 1 
into a tensor product of spinors. This tensor square provides the parameterization 
when restricted to the integer spinors. 

Example. Consider Pythagorean quadruples, that is quadruples of integers (a, 6, c, d) 

a 2 + b 2 + c 2 = d 2 . 

Among the examples are (1,2,2,3), (1,4,8,9), (6,6,7,11), etc. A well-known formula 
that produces Pythagorean quadruples is [10J: 

a — 2mp 

b — 2np 

c = p 2 - (m 2 + n 2 ) ^ 7 ' 2 ^ 
d — p 2j r m 2 + n 2 . 

It is also known that not all quadruples are generated this way, for instance (3, 
36, 8, 37) is excluded [E]. 

Let us try our method. The spinor representation of the Clifford algebra R3 5 i 
is well known. The light-like vectors are represented by the Hermitian 2x2 matrices 
that split into a spinor product, 

d + a b + ci 



M 



ci d 





z 


= 2 




w 



[z w], (7.3) 
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Space of Pythagorean spinors 



E = R 2 (= C) — symplectic space 

(0 = symplectic structure 



A 

n 


i. 




m 



9 



(p([m,n]) = 



Space of triangles 



R 2 1 — Minkowski space 

g = pseudo-Euclidean product 




right 
triangles 

Pythagorean 
triangles 



Symmetries of the structure © 
SlA^R) => SL^Z) => gen { U,R,L 



Symmetries of the structure g 
0(2,R) => 0(2,Z) => gen{U,R,L} 



w 
D 



W 

D 



E 



"(8)" 



-> E®E* <~ 



R 21 



SL^R) 



tt 2 ,2 = spanjai, i=0,..,3} ^~ 
matrix representation 



t 



2,2 



K= R 
Kwaternions 



Figure 7. Objects related to Euclid's parametrization of 
Pythagorean triples. 
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where the matrix M on the left side is well-known in physics with d standing for 
time and a, 6, c for spatial variables. Note that indeed det M = d 2 — a? — b 2 — c 2 
represents the quadratic form of Minkowski space. The factor of 2 is chosen to 
keep things integer. 

Let z = m + ni and w = p + qi. Then (7.3) becomes: 

r ^2 _^^2 

= 2 



d + a b + ci 
b — ci d — a 



ttt + rr 



(mp + nq) + (np — mq)i 



p 2 + q 2 



(mp + nq) — (np — mq)i 
which may be readily resolved for a, 6, c, d. Thus we have proven: 
Theorem 7.1. The following formulae produce all Pythagorean quadruples: 

.2 2 



a 
b 
c 
d 



m 2 + n 2 — p" — q" 
2(mp + nq) 
2(np — mq) 
m 2 + n 2 + p 2 + q 2 



(7.4) 



that is 



(m z + n z — p z — q 2 ) 2 + (2mp + 2nq) z + (2nn — 2mq) z — (m z 



n 



2 i 2\2 

p +q ) 



The quadruple (3, 36, 8, 37) that was not covered by (7.2) may be now 
obtained by choosing (m, q) = (4, 2, 4, 1) 

(4 2 + 2 2 - 4 2 - l 2 ) 2 + (2 • 4 • 4 + 2 • 2 • l) 2 + (2 • 2 • 4 - 2 • 4 • l) 2 = (4 2 + 2 2 + 4 2 + l 2 ) 2 . 

The standard Euclid's parameterization of Pythagorean triples results by 
choosing n = q = (which imposes c = 0), or by m — n and p — q , although this 
time with a redundant doubling in the formulae. Also, choosing only q — will 



result in the system of formulae (7.2). 



Quite similarly, we can treat Pythagorean hexads using the fact that I^i = 
H(2) (quaternionic 2x2 matrices). 

Theorem 7.2. Let us use collective notation m — (mo, mi,m2, m^) andn — (no, ni,n2, n^) 
Also denote in ususal way mn — mono + mini + m2n2 + m^n^. The following for- 
mulae produce Pythagorean hexads: 



(7.5) 



a 


2 i 2 

— m + n 






a\ 


= 2(nomi — 


nim + m 3 n 2 


- m 2 n 3 


a 2 


2(n m 2 - 


n 2 m + min 3 


- m 3 ni 


a 3 


2(n m 3 - 


n 3 m + m 2 ni 


- min 2 


a4 


= 2mn 






a 5 


2 2 

= m — n 







that is for any m,n E t(;e /iave 



2 2 2 2 2 2 

a = a l + a 2 + a 3 + a 4 + a 5 • 
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Proof. A general Hermitian quaternionic matrix can be split into the Kronecker 
product 



ao + as a 


= 2 


P 


a clq — cl§ 




q 



[p 



(7.6) 



where on the left side ao, as E M, and a = a\ + a2^ + a2j + a3&; E H, and on the 
right side q = mo + m\i + m2j + 7713/^ and p = no + n\i + 77,2.7 + 77,3 fc. Resolving this 
equation gives ( |7.5[ ). □ 

For example (ra,n) = ((1,2,2, 1), (2,1,1, 1)) produces 2 2 + 3 2 + 4 2 + 8 2 + 14 2 = 
17 2 . Note that in this parameterization map, tp : Z 8 — > Z 6 , the dimension of the 
parameter space outgrows that of the space of Pythagorean tuples. 



Reconsidering the entries of matrices of equation (7.6) gives us this alternative 
generalization: 



Theorem 7.3. Consider equation 

a + 6 



c 

a — b 



= 2 



P 



[ P 



(7.7) 



with a, b E R, and c,p,q E A, where A is an algebra with a not necessarily positive 
definite norm and cojugation satisfying aa* = |a| 2 and (ab)* = 6* a*. If the qua- 
dratic form of A is of signature (r, s), then formula ( f7.5[ ) produces Pythagorean 
(r + 1, 5 + 1) -tuples. 

Note that as a special case we may use the Cliford algebras themselves as 
algebra A. As an example consider "duplex numbers" D, [7J, which form Clifford 
algebra of R 1 . Let c = Co + c\I E D, where co,ci E R, and I 2 = 1 (pseudo- 
imaginary unit). Similarly, set p — po +P\I and q — qo + Then ( |7.7[ ) produces 
Pythagorean (2,2)-tuples: 

(Po -Pl + Qo- qlf + (2/Wi - 2pig ) 2 = (Po -p\-ql + q\? + (2po^o - 2p iqi f 

A simple application that goes beyond Clifford algebras: using the algebra of oc- 
tonions, A = O, results in a parameterization tp : Z 16 Z 10 of Pythagorean 
"decuples" by 16 parameters. 

Remark 7.4. Since Clifford algebras are 2 n -dimensional, using them in Theorem 



7.3 will lead to parameteriztion of (generalized) Pythagorean (2 n + 2)-tuples by 
T^ 1 "Euclid's parameters". For n = 0, 1, 2, . . . we get 3-, 4-, 6-, 10-, -tuples. 
Incidentally, these numbers occur frequently in various string theories, the reason 
for which is not fully understood [T9] . 
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